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Abstract 

The illuminated dusty surface of Mars acts like a gas pump. It is driven by thermal creep at low pressure within the soil. In the 
top soil layer this gas flow has to be sustained by a pressure gradient. This is equivalent to a lifting force on the dust grains. The 
top layer is therefore under tension which reduces the threshold wind speed for saltation. We carried out laboratory experiments to 
quantify the thickness of this activated layer. We use basalt with an average particle size of 67 yum. We find a depth of the active 
layer of 100 to 200 /im. Scaled to Mars the activation will reduce threshold wind speeds for saltation by about 10%. 
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1. Introduction 

It is a long standing problem how to move particles on the 
martian surface. The most prominent mechanism is wind in 
analogy to transport of sand on Earth. Numerous work has been 
carried out on this in the past especially in wind tunn el experi¬ 
ments dOreelev et al.l.ll98Q[ll99^lwhite et all. 1 1997k . 

Recent images of the HiRISE camera onboard the Mars 
Reconnaissance Orbiter show that martian sand transport is still 
active. They find dunes w hich travel several meters in a few 
years dBridges et al.L l2012h . However, a problem encountered 
in the explanation of particle lift is that wind alone requires 
rather large speed to initiate saltation. The pressure on Mars 
on average is only 6 mbar in contrast to 1000 mbar on Earth. 
This reduces the dynamic pressure of a gas fiow strongly. A 
speed of 30 m/s is supposed to be necess ary to pick up the mos t 
susceptible particles of 100 jim in size dOreelev et al.L Il980h . 
High tangential wind speeds in vortices might also mobilize 
particles. Obviously, dust devils bear witnes s of dust lifting . 
Visible dust devils come in a variety of sizes dkgrenz 20091). 


The la rgest gnes might easily l ift dust. However, IStanzel et al 


d2008h and iReiss et all d20l4 find wind speeds (tangential 
and transversal) w hich are not always large enough. Also, 
Reiss et aP ( 20091) observed dust devil activity on Arsia Mons. 
This relates to an atmospheric pressure of only 2 mbar which 
requires still larger wind speeds. 

There have been suggestions to support or ease particle 
lift one way or the other. The choice of particles to be 
picked up has been varied in wind tunnel experiments. As an 
example, the rolli ng of volcan i c glas s particles might reduce 
threshold speeds dde Vet et all I20l4 . The pressure within 
dust devils is reduced compared to ambient conditions. It has 
been proposed that the traverse of such a pressure minimum 


might be sufficient to lift dust (iBalme & Hagermaml l2006h . 
Last not least and co nnecting to the work presented here, 
IWurm & KraussI ( 2006h found that illumination of a dust bed 
at low pressure pro vides a lift. This was applied to Mars by 
IWurm et ^ (12008 ). Especially this latter effect is strongly 
depending on ambient pressure in a maybe non-intuitive 
way. Wind or gas drag and dynamic pressure decrease with 
decreasing pressure. The induced lifting force of an insolated 
surface can increase to lower pressure by orders of magnitude 
in strength. The force peaks around Knudsen numbers of 
Kn ^ 1, where Kn is the ratio between the mean free path 
of the gas molecules and the size of a particle or pore within 
the dust bed. Hence, for micrometer dust particles insolation 
supported lift is not important on Earth but maximized on Mars. 

The model discussed for this lifting force so far included pho- 
tophoretic forces, solid stat e greenhouse effects and gas com- 
press i gn by thermal creep ( Kocifaj et~^. 1201 ll : iKelling et al.l 
1201 1I: Ide Reiile et al.L 120131) ^ These are important on long 
timescales (hours) as current research is supporting (Koester 
et al. personal communication). However, here we consider 
gas fiow through the dust bed and related pressure differences 
which was not included in the earlier models. This provides 
lift for a sample where illumination changes on short time 
scales of seconds or even fractions of seconds. The importance 
of this became obvi ous in microgravity experiments where 


de Beule et al.1 (120141) observed an efficient gas fiow through an 
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illuminated dust bed directed upwards. The effect is tied to the 
temperature profile within the illuminated dust bed. 

The basic pumping is provided by thermal creep. Eor a con¬ 
stant temperature gradient pure thermal creep pumping does 
not require any pressure differences within the soil. However, 
the temperature profile along the top layer of an illuminated 
soil is rather fiat as radiation is absorbed and thermally emitted. 
This layer cannot pump by thermal creep but has to keep 
up the thermal creep gas fiow from below nevertheless. A 
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pressure gradient is established close to the surface to do so. In 
accordance with Darcy’s law the pressure increase below the 
surface transports the gas flow handed over by the Knudsen 
pump (thermal creep) below. Once set, this overpressure not 
only moves the gas but also acts on the dust particles. Particles 
within this top layer are under constant tension and can be 
ejected if cohesion and gravity can be overcome by any means. 
We call this an insolation activated layer. 

We quantify the thickness of this layer here based on labora¬ 
tory experiments. We illuminate a dust sample with sufficient 
light flux to compensate gravity. We then remove cohesion by 
short vibrations. This leads to a removal of the tension activated 
layer down to the pressure maximum and allows its thickness 
to be determined. 


2. Sub-surface pumping 

We detail the light induced sub-surface pumping in this sec¬ 
tion. The principle is shown in Fig. \T\ 



Figure 1: Principle for pressure distribution for given temperature profile and 
open geometry (ambient pressure the same on both sides of the dust sample). 
Lx 2 marks the depth of the dust bed from 0 to xi with constant temperature and 
Lxi the part of the dust bed with a temperature gradient from xi to xq. 


The light enters the dust bed and is absorbed. The heat 
is conducted further down into the dust bed. In addition, at 
the surface the dust bed can cool by thermal radiation. In 
consequence a temperature gradient is established starting a 
few particle layers within the dust bed (at depth xi in Fig. 
[T]) and is directed to deeper layers. Temperature gradients 
always lead to a thermal creep gas flow in a porous medium, 
where gas is transported along the particl es’ surface from c old 
to warm. This was first introduced by iMaxwe il d 18791) as 


thermal transpiration where two gas reservoirs with different 
temperatures are connected by a small channel in a low 
pressure environment. Gas molecules creep along the channel 
wall from cold to warm. If the diameter of the channel is 
comparable to the mean free path of the gas molecules the 
pressure driven back flow can be smaller than the thermal creep 
flow. In an illuminated dust bed the gas molecules creep from 
cool layers deep within the dust bed (vq in Fig. [T]) upwards 
until the temperature levels off close to the surface at xi . 


If the temperature increase in Fig. \T\ would be linear from 
Xq to the surface there would be no pressure differences. Ev¬ 
ery sub-layer would just transport the same amount of gas by 
thermal creep. However, if there is a top layer of constant tem¬ 
perature the gas molecules don’t creep along these particles ah 
the way to the surface, but only to xi . The thermal creep gas 
flow leads to a concentration of molecules and the pressure is 
locally increased. 

The increase of pressure below the surface at xi leads to a pres¬ 
sure driven gas flow through the top layer. The pressure adjusts 
itself to a value where gas flow through the top layer matches 
the incoming thermal creep gas flow from below. Both aspects, 
the temperature g radient dr i ven th ermal creep gas flow (Knud¬ 
sen compression (lKnudsen.ll909k ) and the pressure driven gas 
flow (Darcy flow dParcvi 1856h ) are usually occuring in the 
same capilla ry. The mass flow of the gas through capil l aries w as 
described bv fSone & Itakural (ll99Qh and iMuntz et al.l (l2002h as 
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where pavg and Tavg are the average pressure and temperature 
within the dust bed, F is a factor giving the amount of capil¬ 
laries within the surface area A, ks is the Boltzmann constant, 
p the molecular mass of the gas, Lr and Lx are the radius and 
length of the capillaries and AT and Ap are the temperature 
and pressure differences within the dust bed, respectively. 
The coefficients Qp and Qp depend on the Knudsen number 
and describe the pressure driven (back) flow and the flow by 
thermal creep, respectively. It has to be noted that the length 
of the capillaries Lx is different in the thermal driven (Lj^, 
Uo ^ -^il) and the pressure driven {Lx 2 , |0 ^ xi|) part. 


If the dust bed is heated two mass flows develop. The first 
one Ml is dominated by thermal creep (Knudsen pump) due to 
the temperature gradient 
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and the second one is driven by the pressure building up accord¬ 
ing to 
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as the top layer temperature is flat or It might 

be noted that the latter is also equivalent to a description by 
Darcy’s law. In total, the flow velocity is set by the Knudsen 
pump within the dust bed and the overpressure which maintains 
that flow also in the top layer activating it by putting tension on 
the dust particles. If the force caused by this pressure gradient 
overcomes gravity and cohesion, particles can be lifted from 
the dust bed’s surface. 


To quantify the temperature profile and the resulting flow ve¬ 
locities and pressure differences we mod eled the insolation o f 
a dust bed as described in earlier work bv iKocifai et al.l ( 2011 ). 
Figure [2| shows the numerical simulation of a temperature pro¬ 
file in a dust bed consisting of 25 yum (radius) spheres. The 
calculations were chosen to closely match the experimental set¬ 
ting: a composition of 25 [im sized (radius) spheres in a dust 
bed w ith a thermal conductiv i ty of 0 .01 W K“^m“^ (at = 10 
mbar ( Presley & Christensenl 1997 )1 and illumination with 4 
kW m“^ (see Kocifaj et~^ (1201 1 ) for details). The temperature 
starts to decrease beyond 100 yum. Therefore, we expect that the 
activated layer would show a similar spatial scale. 



Figure 2: Simulation of temperature profile within a dust bed consisting of 25 
jum spheres. The dust bed has a thermal conductivity of 0.01 W K“^m“f The 
light flux is 4 kW and illumination times are on the order of seconds. 


3. Experiment 

The setup of the experiment is shown in FiglJl Basaltic dust 
with grain sizes ranging from 0 to 125 yum (average size of 67 
yum (diameter) in a volume distribution, see Fig. © is placed 
between two glass plates. It forms a thin dust bed with a width 
of 2 mm and 25 mm height. 
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Figure 3: Schematic setup of the experiment: A dust sample is placed between 
glass plates of 2 mm distance on a table inside of a vacuum chamber with 
ambient pressure of 0.1, 1 and 10 mbar. The dust is illuminated by a red laser 
(655 nm) from above. A weight inside the chamber can be dropped upon the 
table, leading to a tension release within the dust. In addition some particles 
move upon the impact without being illuminated. Images are taken before and 
after the impact within and outside the laser beam. For the images outside the 
laser beam an additional light source is used to observe the surface with the 
camera. 


the dust particles and - if present - releasing the tension due to 
the overpressure. This installation is placed within a vacuum 
chamber and experiments have been carried out at an ambient 
pressure of 0.1 mbar, 1 mbar and 10 mbar. A red laser (655 
nm) illuminates a spot of 8 mm width of the dust sample for ~ 
20 seconds, providing a flux density of 4 kW m“^. This flux 
is chosen to be close to the limit before a continuous particle 
ejection occurs. It provides the most tension to ensure that the 
whole activated layer is ejected on tension release. 


Images of the dust surface before and after the impact are 
taken. After each measurement the laser spot is moved further 
over the surface. After 5 measurements new dust is added to 
the sample. The images within the laser beam show a light 
surface and a black background (Fig. [5l la and lb). In addition 
images outside the laser beam are taken before as well as 
after the impact to measure the particle loss without radiation. 
These are taken as bright held images with transmitted light to 
have a black surface on a bright background (Fig. [Jl 2a and 2b). 


About 0.5 mm of dust emerges above the plates, forming 
a smooth surface which can be observed by a camera and 
microscope optics. 

In a distance of 3 cm to the glass slides a weight of 30 g is 
locked 8 mm above the experimental table by an electromagnet. 
If the magnet is turned off, the weight drops. This results in a 
short vibration of the dust sample, removing cohesion between 


The surface lines before and after an ejection are traced. The 
particle loss is calculated by substracting the lines of the surface 
before and after the impact. For each image an average of the 
thickness of the removed layer with and without illumination is 
determined. A total of 40 to 50 averaged measurements is taken 
for each pressure. 


3 




































Figure 4: Volume distribution of the used sample. 



Ub = 1.38 X 10“^^ J K“^. We assume a temperature differ¬ 
ence Ar = 100 K over a length L^i = 1 mm (see Fig. O 
and an average temperature of Tavg = 273 K - f Ar/2 K within 
the he ated part of the dust bed. According to ISone & Itakural 
( 1990b . the coefficients are Qt = 0.22 and Qp = 1.6. This re¬ 
sults in Apmax = 0-4 mbar. 

In more detail the overpressure ^Pl ^2 equals App^^ in Eq. 0, 
and eventually, the gas velocity is the same in the temperature 
and the pressure driven part, as well as the mass flow. Solving 
the equation Mi = -M 2 for the overpressure Ap, we get 




^jcl + ^jc2 


Pavg 


Ar Qt 

Tavg Qp 


( 6 ) 


Calculating this overpressure we take L^i = 1 mm for the length 
of the capillaries with a temperature gradient and Lx 2 = 0.1 mm 
for the length of constant temperature (see Fig. [2]) we get an 
overpressure of Ap = 0.05 mbar. This is consistently lower that 
the estimate of the maximum pressure difference of Ap = 0.4 
mbar. For an overpressure of about 0.05 mbar, the mass flow 
trough the pores is on the order of 10“^ kg s“^. The velocity of 
the gas flow can be calculated by dividing the mass flow by the 
illuminated area and the density of the gas at the given pressure 
Pgas = P pKTavgRg) with the molar mass of air p = 28.96 AMU 
and the molar gas constant Rg = 8.3 J (K mol)“^. At given pa¬ 
rameters and p = 10 mbar, the gas flow is 11 cm s“^. This is 
consistent with measurements of trac er particles in drop tower 
experiments bv ide Beule et al.l (12014 ). 

Further on, the force induced by the overpressure can be calcu¬ 
lated for a particle column as 


Figure 5: Example of the dust bed within (la) and outside (2a) the laser beam at 
1 mbar ambient pressure. The images lb and 2b show the results after tension 
release. The dust bed within the laser beam (la, lb) was illuminated by a red 
laser (680 nm) with a light intensity of 4 kW m~^. 


4. Model 


The mass flow of gas through capillaries can be calculated 
by Eq. ([T]). We assume that two gas reservoirs are connected 
by a small capillary with a constant temperature gradient. For 
Kn ^ 1 the gas is driven by thermal creep and flows with a 
certain rate through the capillary until the rising pressure dif¬ 
ference Ap balances the thermal creep flow as a pressure driven 
flow through the top layer. There is a maximum Ap aX which 
no net gas flow is active (M = 0, top layer as lid). From Eq. 0 
we get (IMuntz et al. . [2002h 
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We model the dust bed as follows: The mean radius of the 
dust particles is r = 25 pm as used for the temperature profile 
in Fig.0. We consider this a suitable approximation for the 
dust sample of 33,5 pm radius (Fig.0). We further approxi- 
mate the radius of the c apillaries by a pore size = 0.2 x r 
( Jankowski et all l2012h . The ambient pressure is p = 10 


mbar (in one of the experimental sets). The illuminated area 
of the dust bed is A = 2 mm x 8 mm, the Boltzmann constant 


Fp = Ap cr. (7) 

With a particle cross section cr = nr^ (r = 25 pm) we get 
F = 6.1 X 10-9 N. 

The gravitational force which has to be overcome can be ap¬ 
proximated by Fg = M g Lx2(2r)~^ with m = Al?>7ir^ p. With 
the density of basalt p = 2890 kg m"^ we get Fg = 3.1 x lO'^ 
N, which is lower than the force caused by the overpressure 
by a factor of about 2. This is consistent with the fact that 
we adjusted the light flux to a value where no constant parti¬ 
cle ejections occured, i.e. where gravity is compensated but not 
yet cohesion. As the overpressure is sufficient to compensate 
gravity of the top layer, only cohesion need to be overcome. As 
cohesion is removed by the impact the tension release leads to 
the ejection of the whole active particle layer. 

5. Results 

The measurements show a significant difference between the 
depth of particle loss within and outside the laser beam for each 
pressure value as seen in Fig. 0 Fig. 171 Fig. [Hand in Tab. ffl 
This proves that there is an activated dust layer. The green bins 
in the figures show the probability of the depth to which particle 
loss occurs without light. The blue bins show the probability to 
find a certain layer depth within the laser beam. The negative 
depths refer to the fact that the surface can not be prepared per¬ 
fectly smooth and therefore some aggregates can move during 
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Figure 6: Probability to find a certain depth of particle layer ejected at 10 mbar; 
blue bins are within the laser beam; green bins show the depth without radiation 
for comparison 


Figure 8: Probability to find a certain depth of particle layer ejected at 0.1 mbar; 
blue bins are within the laser beam; green bins show the depth without radiation 
for comparison 



depth [pm] 


Figure 7: Probability to find a certain depth of particle layer ejected at 1 mbar; 
blue bins are within the laser beam; green bins show the depth without radiation 
for comparison 


the tension loss, falling onto the measured spot and increasing 
the surface instead of decreasing it. Without light the measured 
values are strongly peaking around a depth within the dust bed 
of 50 yum. Within the illuminated area of the dust bed’s surface 
the measured depth ranges from 50 up to 500 /im for 1 and 10 
mbar, where the variation of the values seems somewhat higher 
for 1 mbar. For 10 mbar the layer size shows less variation and 
is more centered around 100 yum. The values for 0.1 mbar show 
no significant loss attributed to illumination. This is due to the 
decrease in absolute pressure and pressure difference. While 
the temperature profile is similar as in the other two cases it 
is not imprinting itself enough on an overpressure to dominate 
over Earth gravity. 

The temperature profile for the relevant time scales in our ex¬ 
periment shown in Fig. [2] is consistent with the experimental 
results. The fiat temperature region has an extension of around 
100 yum which is the same as the measured depths for 1 and 10 
mbar. Average values are given in Tab. [T] 


pressure 

mean thickness 

mean deviation 

0.1 mbar 

1 mbar 

10 mbar 

33 pm 

166 pm 

109 pm 

+ 32 pm 
± 57 pm 
± 35 pm 


Table 1: The mean thickness of the active dust layers for the different pressures. 
Each data point is based on about 40 measurements and the mean thickness was 
calculated after subtracting the mean offset. 


6. Application to Mars 


The laboratory setting was chosen to quantify the active layer 
depth in a convenient way. On Mars the light flux is and the 
average temperature are lower and the durations of insolation 
are longer. We therefore adjusted our simulation to scale the 
experimental results to martian conditions. We choose 700 W 
m“^ for the solar insolation and illuminate the dust bed in our 
simulation for up to 12.5 hours at an average temperature of T 
= 200 K. The temperature profiles are shown in Fig. [9j The 
maximum temperatures are somewhat larger than temperatures 
found on Mars. This is due to the simplifications of the model 
but we consider the results to capture the essential trend. Addi¬ 
tionally we consider CO 2 instead of air with a molecular mass 
of mco 2 - 44 AMU, a geometric radius of 4.63 x 10“^^ m and 
a resulting cross section a = 1.6 x 10“^^ m^. At p = 6 mbar 
with a particle density n = pk'^^T~^ the mean free path of CO 2 
is T = ( ^/2 ncr)~^ =19 pm Henc e, the ratio between and 
Qp is 1.68 dSone & Itakural (Il990l) ). According to Fig. [9] we 
expect typical temperature differences of AT = 150 K for a 
length Ljci = 50 mm (capillaries with a temperature gradient) 
and Lx 2 = 1 mm for the activated layer of constant temperature. 
Equation [6] gives a pressure difference of Ap = 0.01 mbar. This 
corresponds to a lifting force for a particle pile in the active 
layer of Fp = 1.8 x 10“^ N. Gravity on the same particles is 
FcMars = m 0.38 g Lx2(2r)~^ = 1.4 x 10“^ N. That means that 
the activation accounts for 13% of gravity in a 1 mm thick layer 
of the soil. Variations in strength are expected for changing 
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Figure 9: same as fig. [2|but for a light flux of 700 W m~^, average temperature 
of 200 K. The illumination times are 0.5 h, 1 h, 1.5 h... 12.5 h from bottom to 
top. 

insolation and soil parameters. 

7. Conclusion 

Our laboratory experiments showed that the illumination of 
a dust bed at low ambient pressure leads to an overpressure 
below the surface. This pressure provides a lifting force on 
particles in insolated dust beds. The pressure and thickness of 
this activated layer can well be described by a thermal creep 
model. In the experiments we used 4 kW m"^ and illuminated 
dust beds for about 20 s. This results in an activated layer of 
about 100 micrometer thickness. In the laboratory the over¬ 
pressure at 1 mbar and 10 mbar is large enough to compensate 
(Earth) gravity. 

At martian conditions the activated layer is supposedly deeper 
up to 1 mm. The supporting force can be on the 10% level 
compared to gravity. We expect that this mechanism supports 
other lifting mechanisms significantly, e.g. by decreasing the 
threshold wind speed for saltation by 10%. 
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